Evidence for weak electronic correlations in Fe-pnictides 
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Using x-ray absorption and resonant inelastic x-ray scattering, charge dynamics at and near 
the Fe L edges is investigated in Fe pnictide materials, and contrasted to that measured in other 
Fe compounds. It is shown that the XAS and RIXS spectra for 122 and 1111 Fe pnictides are 
each qualitatively similar to Fe metal. Cluster diagonalization, multiplet, and density- functional 
calculations show that Coulomb correlations are much smaller than in the cuprates, highlighting 
the role of Fe metallicity and strong covalency in these materials. Best agreement with experiment 
is obtained using Hubbard parameters U < 2eV and J ~ 0.8eV. 



PACS numbers: 74.70.Dd,78.70.Dm,71.10.Fd,71.15.Mb 
I. INTRODUCTION 

In the new and rapidly developing field of Fe-pnictide 
superconductivity, the question of what constitutes the 
basic ingredients for high transition temperatures re- 
mains largely unanswered. Parallels have been drawn 
to the cuprate high-temperature superconductors, which 
contain partially filled d-electron spins that in the par- 
ent phase are aligned antiferromagnetically like the pnic- 
tides, and high-temperature superconductivity emerges 
when magnetism can be suppressed. Common to many 
ideas is that superconductivity itself may be emergent 
from the two competing phases, driven by an underlying 
quantum critical point. 

A key question that needs to be addressed to under- 
stand this framework is whether or not the Fe pnictides 
are strongly correlated like the cuprates. Since den- 
sity functional theory (DFT) calculations have indicated 
that the electron-phonon interaction is too weak to ac- 
count for high transition temperatures, ^'^'^ the strength 
of the Coulomb correlations would give some account for 
the pairing strength possible in an electronic-based pair- 
ing mechanism.^ Recent renormalization group flow^ and 
random-phase approximation (RPA) calculations of effec- 
tive tight-binding models fitted to DFT bands^'^ indicate 
that several pairing instabilities, such as sign-changing 
s— wave pairing and type pairing, all have nearly the 
same energy, which depends subtly on Coulomb param- 
eters U and magnetic exchange J. Therefore pinning 
down these numbers would greatly focus the discussion 



of the physics of the Fe pnictides. 

Theoretically, the situation is complicated. Traditional 
DFT methods, which can be extremely accurate in un- 
correlated materials, can account for the correct atomic 
structure but yield large sublattice magnetic moments 
that have not been observed in experiments.^ This over- 
estimation of the magnetic moment is exactly opposite 
to the situation in the cuprates, where DFT underesti- 
mates the moment, and implies that Coulombic effects 
are a small part of the story for the pnictides. However 
it does not provide an explanation as to why the moments 
are so much larger than those found in experiments^. 

Recently, theoretical treatments using combinations of 
DFT and dynamical mean field theory (DMFT) have 
yielded opposite conclusions. ^^'^^'^^'^^'^^'^^ One set of cal- 
culations yield Hubbard parameters [7^4 eV, giving a 
Kondo-like peak near the Fermi level and a well sepa- 
rated lower Hubbard band, and argue that these materi- 
als are on the verge of a Mott transition. ^^'^^ Another set 
however gives U ^ 1 eV.^^'^^'^^'^^ Angle-resolved photoe- 
mission (ARPES) studies have shown strong density of 
states (DOS) near the Fermi level relative to the pnicto- 
gen valence band.^^ The Fe conduction band states, with 
a bandwidth W 4 eV, as well as overall band disper- 
sions can be well matched to density functional calcula- 
tions with a mass renormalization of 2. This is in direct 
contrast with the cuprates, with valence band spectral 
weight spread out over much larger energy scales, and 
conduction bands having much larger band renormaliza- 
tions, smaller bandwidths and DOS. Thus it would ap- 
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pear from ARPES that these materials are categoricahy 
different than the cuprates. 

X-ray measurements have been crucial in uncovering 
the physics of the cuprates, identifying multi-particle 
states, such as the Zhang Rice singlet, as signatures of 
strong level Coulomb interactions, as well as exci- 
tons and satellites due to strong core- hole interactions.^^ 
These strong satellites and spread out spectral weights 
have not been observed in recent x-ray absorption (XAS) 
measurements at the oxygen i^-edge in 1111, setting an 
upper bound on the effective Coulomb parameter U ^ 1 
18,19,20 Absorption and emission studies on the Fe ^2,3 
edges also are in agreement with weak electronic corre- 
lations, and can be simply matched to the unoccupied d 
DOS determined from DFT calculations. However cal- 
culations pertinent to the XAS process, where core holes 
are created, and the emission process, where photons are 
emitted from the valence states in the presence of a core 
hole, have not been carried out. This is crucially needed 
in order to understand the true role of Coulomb interac- 
tion in these materials. 

In this paper we present a comprehensive study of 
XAS measurements and resonant inelastic x-ray scat- 
tering (RIXS) at the Fe ^2,3 edges in a variety of Fe- 
based materials, including the superconducting 1111 Fe 
pnictide SmOo.ssFeAs and the undoped 122 pnictides 
BaFe2As2 and LaFe2P2- It is shown that the XAS spec- 
tra of Fe pnictides look qualitatively, and in some cases 
quantitatively, similar to Fe metal and show no features 
resembling the multiple peak structures seen in Fe insu- 
lators, such as hematite (a-Fe203) and other iron oxides. 
A resonance study of x-ray emission across the L2 and 
I/3 edges demonstrates that the RIXS spectra is domi- 
nated by fluorescence, with no observance of discernable 
excitonic or satellite peaks. 

In addition, we present calculations using three sepa- 
rate models which specifically include and account for the 
role of the core hole in x-ray absorption and emission pro- 
cesses. These calculations are performed using quantum 
cluster, multiplet, and DFT-based methods, and high- 
light the roles of Fe metallicity, FeAs covalency, and local 
Coulomb and Hund's couplings. DFT calculations using 
PEPp2i g[yQ quantitative agreement with XAS measure- 
ments and align absorption peaks to Fe d-DOS above the 
Fermi level, demonstrating the minor role of core hole in- 
teractions. Cluster calculations of XAS support the role 
of strong Fe-As hybridization involving As states below 
the Fermi level, setting an upper bound ofU ^ 2eV. This 
indicates that the FeAs materials are weakly correlated 
and that the physics is governed largely by Fe metallicity. 

The outline of paper is as follows. In Sec. II, XAS and 
RIXS measurements on superconducting SmOo.ssFeAs, 
non-superconducting BaFe2As2, LaFe2P2, Q^-Fe203, and 
Fe metal are presented, comparing and contrasting qual- 
itative behaviors across these compounds, while in Sec. 
Ill calculations are presented for XAS and XES at the 
Fe L2,3 edges. Sees. Ill A and IIIB present calculations 
for L-edge XAS in Fe clusters to highlight the expected 



role of strong Coulomb correlations, and it is shown that 
spectral features related to the strong correlations that 
are not seen in experiments of Sec. II can be used to set 
upper limits on Hubbard parameters. Moreover, DFT- 
based FEFF calculations, which include multiple scatter- 
ing and core hole effects, are presented in Sec. IIIC, 
and are shown to provide excellent agreement with the 
measured XAS spectra. Finally, Sec. IV summarizes our 
findings and states our conclusions. 



II. XAS AND RIXS MEASUREMENTS 

The SmOo.ssFeAs samples with superconducting tran- 
sition temperature (Tc) of 55K, so far the high- 
est Tc in the family of iron arsenides, were pre- 
pared by a high-pressure synthesis method^^'^^. Sam- 
ple quality was checked by x-ray powder diffraction 
and Tc was confirmed by both transport and mag- 
netic measurements^^'^^. We have also measured F- 
doped samples with the same Tc, as well as the 
non-superconducting parent compounds SmOFeAs, but 
found no obvious difference in the spectra. BaFe2As2 
and LaFe2P2 single crystals were prepared by the flux 
method^^'^^'^^. Data shown here were collected at room 
temperature with incident beam 45 degrees to sample 
surfaces. We noticed serious surface oxidization effects 
for the iron pnictides, and to avoid this surface oxidiza- 
tion problem all the data were collected on in-situ cleaved 
sample surfaces. 

XAS and RIXS measurements were performed at 
beamline 8.0 of the Advanced Light Source at Lawrence 
Berkeley National Laboratory. The undulator and spher- 
ical grating monochromator supply a linearly polar- 
ized photon beam with resolving power up to 6000. 
RIXS data were collected by a Rowland circle geometry 
spectrometer^^ perpendicular to the incident beam. The 
linear polarization of the incident beam is parallel to the 
scattering plane. XAS spectra were collected by measur- 
ing sample current (TEY) and fluorescent yield (TFY). 
All XAS spectra have been normalized to the beam flux 
measured by a clean gold mesh. The resolution is better 
than 0.2eV for XAS measurements. For the X-ray emis- 
sion measurements, the incident beam resolution is about 
0.9eV and the spectrometer resolution is about 0.7eV. 

The Fe ^2,3 absorption structure of iron pnictides are 
shown on top of Figs. 1-3. According to dipole selection 
rules, iron is a 3d element displaying L2,3 absorption fea- 
tures from 2p^3d^ to 2p^3d^ transitions. The spin-orbit 
interaction splits the 2p states into 2pi/2 and 2p3/2, lead- 
ing to two well separated peaks. The intensity ratio of 
the two peaks is largely defined by the high-spin or low- 
spin ground states related to the crystal field^^ . As the 2p 
core levels are featureless and narrow, L2 and L3 absorp- 
tion peaks often provide detailed information on the elec- 
tronic structure of the unoccupied 3d states. As shown 
in Figs. la-3a, all the iron pnictide samples, including 
the 55K Tc SmOo.ssFeAs (Fig. la), non-superconducting 
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FIG. 1: (a) X-ray absorption spectra of the 55K Tc 
SmOo.ssFeAs. TEY and TFY intensity of the Fe L2,3 edges 
is plotted as a function of incident photon energy. The dif- 
ference between TEY and TFY data is mainly from the self- 
absorption effect in TFY. (b) RIXS spectra of SmOo.ssFeAs 
collected with excitation energy across the Fe L2,3 absorption 
peaks. The number on the left stands for the excitation en- 
ergy corresponding to the number marked in (a), the value of 
which is marked with the arrows on the spectra. Inset shows 
the prominent Fe L3 emission peak collected with excitation 
energy above the Fe L3 absorption edge (no. 2 to 8), they all 
overlap nicely with the nonresonant spectrum (no. 8). Note 
that RIXS spectra were normalized to the Fe L3 emission 
peak for emphasizing the similar lineshape. 



BaFe2As2 (Fig. 2a) and LaFe2P2 (Fig. 3a), exhibit only 
the two major peaks, L2 at about 720eV and L3 at about 
707eV. There are weak shoulders around 709. 5eV, but 
no peak splitting or intensity ratio change was observed. 
This result is consistent with that on other llll^^and 
122 compounds^^. As XAS has been demonstrated to 
be a powerful tool for probing the crystal field and elec- 
tronic interactions for 3d metals^^'^^'^^; the non-splitting 
XAS structure indicates a weak crystal field effect^^ that 
favors high-spin ground states. 

Fig. lb shows the RIXS data of the superconducting 
SmOo.ssFeAs obtained at energies labeled in Fig. la. The 
spectrum on top (No. 8) was collected with an incident 
photon energy of 735eV, which is far above the Fe L2 and 
L3 absorption edges, the so called nonresonant normal 
emission spectrum. Like all other 3d transition metals, 
this nonresonant spectrum exhibits two main fluorescent 
features at about 704eV and 717eV, resulting from the 
refill of the 2p3/2 and 2pi/2 holes respectively. The 2pi/2 
feature is very weak compared to the L2 edge in the XAS 



FIG. 2: (a) Fe L2,3 XAS spectra of a BaFe2As2 single crys- 
tal, (b) RIXS spectra of BaFe2As2 collected with excitation 
energy labeled and marked in (a). Inset shows that all the 
Fe L3 emission peaks (scaled to the same intensity) collected 
with excitation energy above Fe L3 absorption edge (no. 2 to 
8) overlap with the nonresonant spectrum (no. 8). 



spectrum, partially due to the Coster-Kronig decay pro- 
cess of the 2pi/2 holes to 2^3/2^^- 

The RIXS spectra collected with resonant energies also 
display the strong 704eV peak as seen in the nonreso- 
nant spectrum. With the excitation energy approaching 
the L3 absorption edge (No.l), the 704eV peak evolves 
and stays at the same energy with all the excitations 
above the L3 edge (No. 2-7). This fluorescent feature 
does not track the excitation energy and overlap with 
the 704eV peak in the nonresonant spectrum (inset of 
Fig. lb). No energy loss feature, which is normally asso- 
ciated with various electron excitations and correlations, 
was displayed by the RIXS data. 

RIXS of nonsuperconducting BaFe2As2 (Fig. 2b) and 
LaFe2P2 (Fig. 3b) share the same characterization as 
that of the 55K superconducting SmOo.85FeAs, also 
in agreement with RIXS data reported on another 122 
compound^^. Charge excitation features like Kondo peak 
and lower Hubbard peak are completely absent, and the 
RIXS data is dominated by a peak at 704eV with the 
only difference being the strength of the 701.5eV shoul- 
der. Again, the prominent peaks collected at different 
resonant energies overlap nicely with the fluorescent peak 
in the nonresonant spectrum (insets of Fig. 2b and 3b). 

The absence of excitation induced energy-loss features 
in the RIXS data for all the iron pnictide samples indi- 
cates the weak correlation in this system. It is thus desir- 
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FIG. 3: (a) Fe L2,3 XAS spectra of a LaFe2P2 single crystal, 
(b) RIXS spectra of BaFe2P2 collected with excitation energy 
labeled and marked in (a). Inset shows that all the Fe L3 
emission peaks (scaled to the same intensity) collected with 
excitation energy above Fe L3 absorption edge (no. 2 to 7) 
overlap with the nonresonant spectrum (no. 7). 
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FIG. 4: (a) Fe L2,3 XAS spectra of Fe metal, (b) RIXS spec- 
tra of Fe metal collected with excitation energy labeled and 
marked in (a). Inset shows that, with the elastic peak track- 
ing excitation energy, all the Fe L3 emission peaks collected 
with excitation energy above Fe L3 absorption edge (no. 2 to 
7) overlap with the nonresonant spectrum (no. 7), same as 
that of iron pnictides. 



able to compare the iron pnictides with known metallic 
and insulating iron components, to reveal the importance 
of metallicity and to show that this result is not exper- 
imental resolution limited. In Fig. 4, we show the XAS 
and RIXS data collected on pure iron metal. The XAS 
(Fig. 4a) displays non-splitting L2 and L3 peaks^^'^^'^^ 
with a slightly weaker shoulder compared to the iron 
pnictides. The RIXS data of Fe metal show more sym- 
metric peaks without shoulders, as well as stronger elastic 
peaks tracking the excitation energies. But just like the 
iron pnictides, the RIXS lineshape is dominated by the 
peak at 704eV which overlaps with the fluorescence peak 
collected with off resonance excitation energy (inset of 
Fig. 4b), and the iron metal resembles all the iron pnic- 
tides in the featureless RIXS data without excitation or 
correlation peaks. 

On the contrary, the a-Fe203 powder sample dis- 
plays very different XAS and abundant features in RIXS. 
The XAS data (Fig. 5a) shows strong splitting struc- 
ture on both L2 and L3 absorption edges due to the 
interplay of crystal-field (10I^g=0.88eV) and electronic 
interactions^^. RIXS data (Fig. 5b) show that the spec- 
tral appearance changes drastically with excitation ener- 
gies, and obviously do not overlap with the nonresonant 
spectrum (No. 10). We plotted the energy loss features 
at different resonant energes in Fig. 5c. Particular en- 
ergy loss features, as indicated by the gray lines, were 



enhanced at particular resonant energies, leading to very 
different lineshape. These energy loss features are sig- 
natures of (i(i-excitations^^, details on which is not the 
topic of this paper. With better resolution, our RIXS 
data revealed more excitation modes than that in previ- 
ous publications, covering the whole range of the optic 
absorption bands^^. 

For a direct comparison between all the samples, and 
for comparing with the theoretical calculations. Fig. 6 
shows the RIXS data collected with 708eV excitation en- 
ergy as well as the XAS data on the five different samples. 
There are only minors difference in the symmetry of the 
lineshape and strength of the shoulders between the XAS 
data of iron pnictides and iron metal; while the crystal 
field splitting leads to very different XAS spectrum of a- 
Fe203 data (Fig.6b). The RIXS (Fig.6a) of iron pnictides 
and iron metal is dominated by the prominent fluores- 
cent peak with no energy loss feature related to charge 
excitations; while for a-Fe203, RIXS shows strong en- 
ergy dependence and complex energy loss structure from 
electronic excitations and correlations. The similarity on 
the spectra between iron pnictides and iron metal, as 
well as the absence of charge excitation features in the 
RIXS data, suggests that iron pnictides are unlikely to be 
strongly correlated systems, which is further elaborated 
by theoretical calculations below. 
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III. CALCULATIONS 

In order to understand the main features of the ex- 
perimental spectra, namely that the spectra of the Fe 
pnictides greatly resemble that of Fe metal, we proceed 
in three steps. First, to determine the importance of cor- 
relations, we present exact diagonalization calculations of 
a Hubbard model cluster which can be solved for either 
strong or weak Coulomb correlations. These calculations 
are used to give a qualitative estimate of the size of the 
Hubbard U and Hund's J. Knowing the values of these 
parameters relative to the band width allows us to deter- 
mine whether weakly-correlated methods are appropriate 
for describing the experimental spectra. Second, to bet- 
ter understand the relationship between multiplet, spin- 
orbit, and crystal field effects on the spectra we present 
atomic multiplet calculations. And finally, having estab- 
lished the relatively minor role of correlations in the Fe 
pnictides, we present DFT-based calculations of XAS and 
XES spectra for comparison with experiment. 



FIG. 5: (a) FeL2,3 XAS spectra of a-Fe2 03 powder, (b) RIXS 
spectra of a-Fe203 collected with excitation energy labeled 
and marked in (a), (c) Energy loss features corresponding to 
(ic/-excitations marked with gray lines. 
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FIG. 6: (a) Comparison of RIXS data on the noted samples 
at 708eV excitation energy. All iron pnictides and iron metal 
show only fluorescent peaks same as the nonresonant spec- 
trum (insets of Fig.l-4b), while a-Fe203 displays multiple 
peaks as the signature of c?(i-excitations. (b) Comparison of 
XAS data on the same five samples, a-Fe2 03 again displays 
very different lineshape due to the crystal field splitting. 



A. Cluster Diagonalization 

In order to see explicitly how correlations affect the 
XAS profile, we perform a model many-body calculation 
based on the exact diagonalization (ED) technique, with 
a multi-orbital Hubbard model as the effective Hamil- 
tonian. This approach has been successfully applied to 
understand the correlated physics in materials such as 
the high-Tc cuprate parent compounds^^'^^'^^. 

The Fe pnictides have a tetrahedral FeAs4 plaquette 
serving as the building block of the two dimensional 
Fe2As2 layer. We have therefore attempted to capture 
the essential physics revealed from XAS spectra with a 
FeAs4 tetrahedral cluster including the five Fe 3d levels 
and the As ^Px,y,z orbitals. The energy eigenstates that 
are necessary for calculating the XAS cross-sections via 
Fermi's golden rule are then obtained by diagonalizing 
the multi-orbital Hamiltonian. 

Our cluster calculations have been carried out in an 
assumed high-spin state for the Fe 3(i-levels, which 
is energetically preferred over the low-spin configuration 
due to Hund's coupling. While the experimentally mea- 
sured magnetic moment in Fe pnictides is about 0.35/iB, 
LDA predicts a larger magnetic moment^. This strength 
would decrease if the system is more delocalized. 

The multi-orbital Hamiltonian entering the calcula- 
tions can be written as = Hk -\- -\- He + Hq . Here 
Hk is the kinetic energy term: 

Hk = ^ tpd,^{d\^PjYa ^ h.c.) 

where dt^^ creates a particle with spin a in orbital 7 at 
the Fe site, and P^^^cr creates a particle with spin a in 
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orbital 7' at As site j. The relations among the multi- 
orbital hoppings are derived from the the Slater-Koster 
table"^^, with the strengths of these Slater-Koster ma- 
trix elements being \Vpda\ =14.091 (eV- A) ^^jj"^ , and 
iVpdTrl = -^I'^dcrl- Here d is the Fe-As bond length (in 

units of A), and the material specific values and 
are either given or can be calculated from Ref.^^ . In this 
work we use the values: d = 2.39A, = 13. 2A, and rd = 
O.744A. We have further assumed that \Vppc 
is the orbital site-energy term: 



2 l^pdo 



He = ^ed{j)nd. 



70" 



(2) 



with Ud 



,7(7 



— ^7cr^7cr7 and Tip j^fj — Pj-ycrPj'jcr' The Fe 



Cg and t2g orbital site energies are defined with respect 
to their center of gravity Cd by €d{eg) = Cd — QDq^ and 
ed{t2g) = + 4:Dq. The arsenic p orbital site energy 
Cp is defined by A = Cd — €p -\- nil for the d^ con- 
figuration, where A is the charge transfer gap energy. 
The subtraction of an average Coulomb repulsion term, 
U = A — -\- ^C^ ensures a d'^ ground state; A, 5, 
and C are the Racah parameters. 



The correlated physics is introduced directly from the 
Coulomb interaction term, including intra-orbital on-site 
Coulomb interactions, Hund's exchange coupling, and 
electron pair hopping processes, written as^^: 



= — ^ nd,-fCTnd,-fCT' + — ^ rid,^and,Ya' + - ^ d^^d^^^^^d^a'dya + y ^ d\^^d^^^,d^'a'd^'cT'{^) 



The above tight-binding parameters are related via the 
Goodenough-Kanamori- Anderson relation: U = /7' + 2J, 
and J = J' . Written in terms of the Racah parame- 
ters, the on-site intra-orbital Coulomb repulsion U is ex- 
pressed as /7 = A + AB + 3C On the other hand, the 
Hund's coupling J, typically of order ~ 1 eV in cuprates 
and other transition metal complexes, are solely deter- 
mined by B and C. Later we shall treat A as an adjust- 
ing parameter, and hence the on-site Coulomb repulsions, 
seeing how it affects the XAS spectra. 

For the XAS final states we include an additional core- 
hole potential term in the Hamiltonian: 

Hq^^UqUd^^an'', (4) 

7,cr 

where = d\dc^ and d\ is the creation operator for a 
hole in the core- hole orbital. The strength oiUq can be 
determined experimentally from the energy separation of 
the well-/poorly-screened resonances, and is of the same 
order of magnitude as U. Here we use \Uq\ = U for 
simplicity. In short, the tight-binding parameters used 
in the calculations are (in units of eV): Vpdo- = —1.10, 
VpdT. = 0.63, Vppcr = 0.55, and Vpp^ = -0.15; B = 0.10, 
and C = 0.40 (resulting in a J {eg) = 0.8); = 0.00, 
10L>g=0.20, and A = 1.50. 

Fig. 7 shows the Fe-pnictide La-edge XAS spectra from 
the cluster calculation, with varied Coulomb repulsion U . 
A stronger U suppresses the XAS shoulder peak inten- 
sity. The shoulder peak is further split into a two-peak 
structure for larger Hund's coupling, as is shown in Fig. 
8. By either increasing the covalency or reducing the 
correlation effects a more featureless XAS spectrum with 
a shoulder peak intensity comparable to experiments is 
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FIG. 7: Fe-pnictide Ls-edge XAS spectra obtained from small 
cluster diagonalization for a fixed J{eg) = 0.8 eV and various 
U. A strong Coulomb repulsion tends to suppress the XAS 
shoulder peak intensity. The dashed line sketches the energy 
separation of the dominant peak and its shoulder . The inset 
is a plot for the peak-shoulder energy separation versus the 
on-site repulsion [/, from which a naive upper bound of [/ ^ 2 
eV is drawn. 



obtained. A naive upper bound for the on-site repulsion 
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FIG. 8: Fe-pnictide Ls-edge XAS spectra obtained from small 
cluster diagonalization for a fixed U = S.6 eV and two differ- 
ent values of J{eg). The XAS shoulder peak is split into a 
two-peak structure by a larger J, indicated by the black ar- 
rows. 
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FIG. 9: X-ray absorption spectrum: Iron L-edge multiplet 
calculation for BaFe2As2. 



is therefore drawn by looking at the XAS shoulder peak 
structure, as well as its energy separation from the domi- 
nant peak. According to the calculation, we estimate the 
Coulomb interactions to he U ~ 2 eV and J {eg) = 0.8 
eV, consistent with the empirical formula^^. This result 
suggests that it is more appropriate to treat Fe-pnictides 
as weakly-correlated systems. 

A limitation of the cluster approach is that while it in- 
cludes interaction between the 3d electrons explicitly, the 
states obtained from the cluster diagonalization do not 
include the atomic multiplet structures associated with 
the spin-orbit coupling of the Fe 2p core- hole. To test 
how this multiplet structure affects the resultant spec- 
tra, we have also computed XAS profiles including the 
atomic multiplets, at the expense of removing the pnic- 
tide atoms from the cluster. 



B. Multiplet Calculation 

X-ray absorption spectra are calculated using Fermi's 
golden-rule, with a finite lifetime for the core- hole. Thus, 
the x-ray absorption intensity may be written explicitly 
in terms of a sum over states Itjji) as 

/H = El(^-Wo)P (5) 

i ^ ^ 

We considered the specific case of XAS experiments 
on BaFe2As2. During the optical transitions, the states 
l^/^o), IV^i) and \ipf) belong respectively to the configura- 
tions 2/3d^, 2p^3d^ and {2p^3d^y. For each of these 
configurations, we have to consider electron-electron in- 
teractions, spin-orbit and crystal-field on an equal footing 
while the radial wave-functions are determined by solv- 
ing the Dirac equation. This leads to a splitting of the 
shells into multiplet levels. To be specific, the crystal- 
field is expressed as an electronic potential of external 
point-charge ions interacting with the considered Fe-ion. 



From there we can evaluate the associated XAS spectra 
within the dipole approximation. 

The obtained XAS spectrum of Fig. 9 shows good 
agreement with the experimental data and mostly ex- 
hibits the L2-L3 splitting coming from the spin-orbit 
splitting of the 2p core- levels. However, the present 
calculation does not involve charge fluctuations, the ex- 
plicit inclusion of these would give rise to satellite peaks 
in a RIXS spectrum and the inclusion of the ligands 
would give rise to the formation of bands where electron- 
electron scattering would occur. 

In the limit of a strongly correlated regime the initial 
and final states cannot be described by a single-site ap- 
proach, hence the XAS spectra would exhibit additional 
peaks that cannot be captured by the present multiplet 
picture. The experimental data do not exhibit such a 
signature: besides the spin-orbit splitting the XAS spec- 
trum is rather featureless and the additional shoulders 
can be explained by the effects of the crystal-field on the 
atomic- multiplet. These facts combine to suggest that 
BaFe2As2 is rather not a strongly correlated material. 



C. DFT-based FEFF Calculation 

For weakly correlated materials XAS and XES spec- 
tra can reliably be computed using ab initio methods. 
Such calculations, however, have a considerable degree of 
complexity because in XAS/XES a core- hole is present in 
the final/initial state. The computer code FEFF^^ is well- 
known for treating the core-hole potential with a high 
level of accuracy. Here we have used FEFF to calculate 
the XAS near the L2 and L3 iron edges, the angular- 
momentum projected density of states (LDOS), and the 
XES spectra for the L3 edge. 

Our calculations begin by overlapping relativis- 
tic Dirac-Fock atomic potentials via the Mattheiss 
prescription^^'^^. This prescription fixes a "Norman" 
radius about each atomic site which contains Zi elec- 
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FIG. 10: The Fe L3 edge XES signal (a) calculated using 
FEFF for the four different metallic iron-containing materials 
considered in the experimental section. The Fe L2,3 edge XAS 
(b) calculated using FEFF for the four different iron- containing 
materials considered in the experimental section. 



trons, where Zi is the atomic number of the atom at site 
i. The overlapped Mattheiss potentials are then used 
as the starting point of a self-consistent (SCF) potential 
calculation which uses the ground-state von Barth-Hedin 
exchange-correlation potential^^ on all iterations. Given 
the SCF potential, the relativistic radial wave-functions 
and phase shifts associated with each atomic scattering 
site can be calculated. The single-electron Green's func- 
tion for the entire system may then be written using a ba- 
sis of these radial wavefunctions and spherical harmonics 
with system (cluster) dependent coefficients that are cal- 
culated within multiple scattering theory.^^ Finally, the 
XAS can be calculated from Eq. (5), using the Green's 
function to implicitly sum over states: 

/i - -lm{^o\d^G{Ec + nLo)d\^o) , (6) 

where d is the single-electron dipole operator, G is the 
photoelectron Green's function, and the state \ipo) is the 
core state of interest which, in this work, is either the L2 
or L3 edge of iron. The energy of the absorbed photon 
is hco^ and the ~ symbol means that we have neglected 
to write a number of constant prefactors as well as a 
broadened step function limiting the XAS spectrum to 
uj > \Ec\ + fi, where /i is the Fermi level. The XES can 
be calculated from a similar formula, but is limited by the 
complementary step function to the absorption case. The 
presence of the core-hole on the absorbing atom, as well 
as the effect of the Hedin-Lundqvist^^ self-energy (in the 
"plasmon pole"^^ approximation), are also included in 
the FEFF calculations.^^ The LDOS for each type of atom 
is calculated by integrating the spacially- and energy- 
dependent density about each atom within the Norman 
sphere. This allows for an unambiguous definition of the 
LDOS for each type of atom. 

In the panel b) of Fig. 10 we show our FEFF calcu- 
lations of the iron ^2,3 edge XAS for the four metallic 



FIG. 11: The local angular- momentum projected densities 
of states for each of the metallic materials considered in the 
experimental section. Only the largest angular-momentum 
contributions from each type of atom are shown; the p-DOS 
is shown for As, O, and P; the d-DOS is shown for Ba, La, 
Sm, and Fe. We note that the Fe DOS shown is that of the 
absorbing atom with the core hole. 

iron-containing materials considered in the experimen- 
tal section. In the panel a) of Fig. 10 we show our 
XES calculations at the iron L3 edge. These XES cal- 
culations may be compared to the RIXS data for high 
fixed incoming energy and detected photon energies be- 
low \El3 I + /i. Similarly to our interpretation of the XAS 
as reflecting the unoccupied dDOS, the near edge XES 
may be simply interpreted as a reflection of the occupied 
dDOS which is dominated by the Fe contribution near 
the Fermi level. The similarity of the pnictide spectra to 
that of ordinary Fe metal underscores the importance of 
metalicity in these materials. The calculated spectra in 
both the pnictides and iron metal mimic the DFT den- 
sity of states as expected in a weakly-correlated picture; 
the XAS (which involves an initial p-state) is determined 
in the near edge region largely by the unoccupied dDOS. 
The dDOS, in turn, is dominated near the Fermi level 
by the contribution from iron for all the materials con- 
sidered, as shown in Fig. 11. Other features of both the 
XAS and the XES can be matched to the peaks in the 
dDOS shown in Fig. 11. For completeness we also show 
the pDOS of As, O, and P, in Fig. 11, but we note that 
it is not directly related to the XAS or XES spectra. 

The fact that in FEFF the features of the XAS and the 
XES spectra can be matched to the peaks in the Fe 3d 
densities of states suggests that the core-hole is relatively 
unimportant and its potential weak. We can test this 
conclusion by computing the same spectra with a method 
that disregards the potential of the 2p core-hole, but re- 
tains the correct band-structure dipole matrix elements. 
For this we used the plane-wave based DFT computer 
program WIEN2k^^ The resulting XES and XAS spec- 
tra for three types of iron based materials are shown in 
Fig. 12. The splitting and relative intensity cannot 

be determined with this method and are therefore intro- 
duced by hand. The computed XAS and XES spectra of 
the pnictides agree well with both our FEFF calculations 
and the experiments. We observed that in both FEFF and 
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FIG. 12: The iron L2,3 edge XAS and XES, calculated using 
WIEN2k, for three different iron containing materials. 

WIEN2k calculations, the effect of oxygen doping and also 
the effect of replacement of one rare earth with another is 
rather small since the XAS is determined largely by the 
local environment of the absorber which is dominated by 
iron. We also note that the theoretical spectra presented 
here are polarization averaged; we have investigated po- 
larization dependence in our FEFF calculations but find 
no significant differences between polarized and averaged 
spectra in the pnictides. 

IV. SUMMARY 

In the young field of iron-pnictide superconductors 
there are currently several open issues. For instance, it 
is not well understood why the observed ordered mag- 
netic moment is so small. Also the role of local Coulomb 
interactions is not well characterized, as is the degree of 
As hybridization with Fe orbitals near the Fermi level. 
The tendency towards the formation of large iron mo- 
ments due to local Hund's rule exchange and a possible 
emerging role for orbital degrees of freedom are other dis- 
puted issues^^'^^'^^'^^. These disputes stand in the way 
of a consensus on the minimal model needed to describe 
the physics and ultimately the pairing mechanism in Fe- 
pnictide superconductors. In this context we have inves- 
tigated five iron containing materials including 122 and 
1111 Fe-pnictides with a combination of XAS and RIXS 
techniques. The first general observation is that the ex- 
perimental data for the Fe pnictides is qualitatively simi- 
lar to other metallic Fe materials and significantly differ- 
ent from large gap Fe-based insulators such as hematite 
a-Fe2 03. 

The three main theoretical approaches that we have 
used to analyze the data incorporate electronic correla- 
tion effects due to electron-electron interactions to dif- 
ferent degrees. If the pnictides were very localized one 
would expect that the essential features of XAS and RIXS 
spectra can readily be captured in a small FeAs cluster. 
Our exact diagonalization computation of the absorption 



spectra for such a cluster in the localized, strong coupling 
limit (large Hubbard U) clearly shows the appearance of 
a high energy peak well separated from the main ab- 
sorption line at the L3 edge. In the experimental data 
this peak is absent - or rather appears as a shoulder of 
the main XAS line. From a comparison of the energy 
position of this shoulder in the data and the cluster sim- 
ulation we extract an upper limit of the Hubbard U of 
2 eV -substantially smaller than the Fe 3d bandwidth. 
Hund's rule J is about 0.8 eV. Coulomb correlations are 
thus much weaker then in the cuprates. This result is 
confirmed by our multiplet calculation, which is a local 
approach that has the advantage of including the spin- 
orbit and Coulomb interactions related to the core-hole. 
The calculated L3-L2 edge energy splittings and intensity 
ratios agree with the experimental data. 

The inference that the Hubbard U is small and the 3d 
electrons weakly correlated, suggests a comparison of the 
spectroscopic data with the results of single-particle, ab 
initio calculations. Such approaches are complicated by 
the fact that a core-hole is present in the final state of 
XAS and the intermediate state of RIXS. The Coulomb 
interaction of the core- hole with the valence electrons can 
in principle be strong and such an electronic correlation 
effect has a profound influence on XAS and RIXS spec- 
tra. The density-functional based FEFF code treats the 
effects of the core-hole potential with a high level of ac- 
curacy. The spectra computed with FEFF confirm the 
presence of the XAS shoulder and agree with the experi- 
mental XAS and RIXS data very well. We have used this 
agreement to further filter out the core-hole induced cor- 
relations. When calculating the spectra with plane-wave 
based WIEN2k code -which includes the proper dipole 
transition matrix elements, but lacks the final state XAS 
core-hole potential- we observe that the simulated spec- 
tra basically do not change. We conclude that in the 
Fe pnictides not only the Hubbard U but also the core- 
hole potential is therefore heavily screened. The elec- 
tronic correlations that the core-hole induces are thus 
weak and consequently the spectra can safely be inter- 
preted in terms of single-particle densities of states and 
the appropriate dipole transition matrix elements. The 
present spectroscopic data and its theoretical description 
thus emphasize the role of strong covalency and Fe met- 
alicity in the Fe pnictides. 
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